Oscillations of complex networks 
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A complex network processing information or physical flows is usually characterized by a number of macro- 
scopic quantities such as the diameter and the betweenness centrality. An issue of significant theoretical and 
practical interest is how such a network responds to sudden changes caused by attacks or disturbances. By 
introducing a model to address this issue, we find that, for a finite-capacity network, perturbations can cause 
the network to oscillate persistently in the sense that the characterizing quantities vary periodically or randomly 
with time. We provide a theoretical estimate of the critical capacity-parameter value for the onset of the network 
oscillation. The finding is expected to have broad implications as it suggests that complex networks may be 
structurally highly dynamic. 

PACS numbers: 89.75.-k, 89.20.Hh, 05.10.-a 
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The response of a complex network to sudden changes such 
as intentional attacks, random failures, or abnormal load in- 
crease, has been of great interest 

since the discoveries of the small-world 1 1 1 ] and the scale-free 
lH2ll topologies. The issue is particularly relevant for scale- 
free networks that are characterized by a power-law degree 
distribution. For such a network, generically there exists a 
small set of nodes with degrees significantly higher than those 
of the rest of the nodes. A scale-free network is thus robust 
against random failures, but it is vulnerable to intentional at- 
tacks 1 1]. This is particularly so when dynamics on the net- 
work is taken into account, which can lead to catastrophic 
breakdown of the network via the cascading process ^ H| 
even when the attack is on a single node. A basic assumption 
underlying the phenomenon of cascading breakdown is that a 
node fails if the load exceeds its capacity. As a result, the load 
of the failed node has to be transferred to other nodes, which 
causes more nodes to fail, and so on, leading to a cascade of 
failures that can eventually disintegrate the network. 

There are situations in complex networks where overload 
does not necessarily lead to failures. For instance, in the In- 
ternet, when the number of information-carrying packets ar- 
riving at a node exceeds what it can handle, traffic congestion 
occurs. That is, overload of a node can lead to the waiting 
of packets but not to the failure of the node. As a result of 
the congestion, traffic detour becomes necessary in the sense 
that any optimal routes for new packets on the network try 
to avoid the congested nodes. This is equivalent to a change 
in the "weights" (to be defined more precisely below) of the 
congested nodes and, consequently, to changes in the macro- 
scopic characterizing quantities of the network. This situation 
usually does not occur when the network is in a normal oper- 
ational state, but it becomes likely when sudden disturbances, 
such as an attack or an abrupt large load increase, occur. A 
question is then whether the network can recover after a finite 
amount of time, in the sense that its characterizing quantities 
restore to their original values. 

In this Letter, we study a class of weighted scale-free net- 



works, incorporating a feasible traffic-flow protocol, to ad- 
dress the above question. In the absence of any perturbations, 
the network is assumed to operate in its "normal" state so that 
its macroscopic characterizing quantities are constants. We 
find that, after a large perturbation, the network can indeed 
recover but only for large node capacities. When the node ca- 
pacities are not significantly higher than their loads in the nor- 
mal state, a surprising phenomenon arises: the macroscopic 
quantities of the network are never able to return to their un- 
perturbed values but instead, they exhibit persistent oscilla- 
tions. In this sense we say the network oscillates. More re- 
markably, as the node capacities are decreased, both periodic 
and random oscillations can occur. The striking feature is that 
the oscillations, periodic or random, are caused solely by the 
interplay between the network topology and the traffic-flow 
protocol, regardless of the network parameters such as the de- 
gree distribution and the overall load fluctuations. For fixed 
network parameters, the oscillations exist despite of the ex- 
plicit form of the local node dynamics, given the simple rule 
that it "holds" and causes the traffic to wait when overloaded. 
This may have wide implications to many traffic problems. 
For instance, it can provide an alternative explanation, from 
the dynamical point of view, for the recently observed ran- 
dom oscillations in real Internet traffic flow ltl3l[l4ll and give 
some insights to the self-similar oscillations of the traffic flux 
observed in WWW lIllfTrjl. 

We begin by constructing a scale-free network of N nodes 
using the standard growth and preferential-attachment mech- 
anism [12]. We next define the node capacity by using the 
model in Ref. H, 



d = (H-a)L<(Q), 



(1) 



where i, (0) is the initial load on node i, which is approxi- 
mately the load in a normal operational state (free of traffic 
congestion), and a > is the capacity parameter. The load 
Li can be conveniently chosen to be the betweenness IItI il8f1 . 
which is the total number of optimal paths |19] between all 
pairs of nodes passing through node i. To define an optimal 
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path at time t, say at this time the weights associated with 
node i and with node j are Wi(t) and Wj[t) (to be defined 
below according to the degree of traffic congestion), respec- 
tively, where there is a direct link /y between the two nodes. 



The weight of the link is then dij(t) = [wi(t) 



/(*)]/2- 



Given a pair of nodes, one packet generating and another re- 
ceiving, the optimal path is the one which minimizes the sum 
of all weights d%j of links that constitute the path. Finally, 
we define a traffic protocol on the network by assuming that, 
at each time step, one packet is to be communicated between 
any pair of nodes. There are thus N(N — l)/2 packets to be 
transported across the whole network at any time. When a 
packet is generated, its destination and the optimal path that 
the packet is going to travel toward it are determined. 

In a computer or a communication network, a meaningful 
quantity to characterize a link is the time required to transfer 
a data packet through this link. When the traffic flow is free, 
it takes one time unit for a node to transport a packet. When 
congestion occurs, it may take a substantially long time for a 
packet to pass through a node and hence any link from this 
node. For instance, suppose at time t there are Jj (t) packets 
at node i, where Ji(t) > Ci. Since the node can process Ci 
packets at any time, the waiting time for a packet at the end of 
the queue is l+iat[Ji(t)/Ci], where int[-] is the integer part of 
the fraction in the square bracket. These considerations lead 
to the following definition of instantaneous weight for node i: 



Wi(t) 



int 



for 



,N, 



(2) 



from which the instantaneous weights for any link in the net- 
work and hence a set of instantaneous optimal paths can be 
calculated accordingly. For free traffic flow on the network, 
we have J,; (t) < C,; and hence tOj (t) = 1 so that the network 
is non-weighted. In this case, the optimal path reduces to the 
shortest path. 

The above model of traffic dynamics on a weighted net- 
work allows us to investigate the response of the network to 
perturbations in a systematic way. In particular, since the node 
capacities are the key to the occurrence of traffic congestion, 
it is meaningful to choose the capacity parameter a in Eq. 
U) as a bifurcation parameter. To apply perturbation, we lo- 
cate the node with the largest betweenness B max in the net- 
work and generate a large number of packets, say ten times of 
B m ax, at time t = 0. The network is then allowed to relax 
according to our model. Initially, because of the congestion 
at the largest-betweenness node caused by the perturbation, 
its weight assumes a large value. As a result, there is a high 
probability that the optimal paths originally passing through 
this node change routes. This can lead to a sudden increase 
in the network diameter, which is the average of all optimal 
paths. As time goes by, the congestion will cascade to other 
nodes which adopt the detoured optimal paths and, as a re- 
sult, the diameter will increase very quickly and reaches its 
maximum at some instant. During this process the congestion 
situation is released at the attacking node while it get worse at 
nodes that paths detour to. After this, the network begins to 
"absorb" the congestions according to the load tolerance, i.e. 
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FIG. 1: For a scale-free network of 1000 nodes and for a = 0.4, 
time evolutions of three macroscopic quantities: (a) the normalized 
network diameter, (b) the normalized betweenness centrality, and (c) 
the number of jammed nodes. The quantities return to their respec- 
tive steady-state values after a brief transient. 
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FIG. 2: (a-c) For the same scale-free network but for a — 0.31, 
periodic oscillations of the normalized diameter (a), the normalized 
betweenness centrality (b), and the number of jammed nodes (c). (d- 
f) random oscillations of the same set of quantities for a = 0.2. 



the capacity parameter, and the recovery process starts. Dur- 
ing this process the network congestion is gradually released 
and the diameter is expected to be decreasing. The same pro- 
cesses apply to other macroscopic characterizing quantities of 
the network, such as the betweenness centrality. However, 
due to the imbalance of load distribution and the high density 
of optimal paths, the final state where the system recovers to 
is highly unpredictable |23]. Therefore we are interested in 
whether these quantities can return to their "normal", or the 
steady-state values before the perturbation. 

For relatively large value of a, the ability of the network 
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guaranteed if aLj(0) > ALj but traffic congestion occurs 
if aLj (0) < ALj. The critical value a c is then given by 
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FIG. 3: (Color online) For a scale-free network of N = 3000 and 
(k) — 4, algebraic scaling of the initial load with the degree vari- 
able k (upper data set) and the scaling of the load change caused 
by unit weight change (lower data set). The parallelism of the two 
sets validates the use of Eq. J3j- The insets show similar plots but 
for different network parameters. These results suggest that the crit- 
ical value a c for network oscillation is independent of the structural 
details of the network. The results are averaged over 50 realizations. 



to process and transport packets is strong, so we expect the 
network to be able to relax to its unperturbed state. This is ex- 
emplified in Figs. [Da-c), the time evolutions of three macro- 
scopic quantities, the normalized diameter (D), the normal- 
ized betweenness centrality (B), and the number of jammed 
nodes nj, respectively, of a scale-free network of 1000 nodes 
for a = 0.4. [For this network the values of the diameter 
and of the betweenness centrality in the unperturbed state are 
(D ) ~ 5.18 and (Bo) w 2.35 x 10 6 . The plotted quanti- 
ties in (a) and (b) are normalized with respect to these "static" 
values.] We see that, after about 7 time steps, these quantities 
reach their maximum values and, after another about 5 steps, 
these quantities return to their respective unperturbed values. 
In this case, the large perturbation causes the network to os- 
cillate but only for a transient time period. As the capacity pa- 
rameter a is reduced, a remarkable phenomenon occurs: after 
an initial transient the network never returns to its steady-state 
but instead, it exhibits persistent oscillations. Figures |2ja-c) 
show periodic oscillations for a = 0.31, where the legends 
are the same as for Figs. Qa-c), respectively. The oscilla- 
tions are in fact period-2 in that each macroscopic quantity 
can assume two distinct values, neither being the steady-state 
value, and the quantity alternates between the two values. For 
smaller value of a, random oscillations [20] occur, as shown 
in Figs. Etd-f) for a = 0.2. 

The critical value a c of the capacity parameter, below 
which persistent network oscillations can occur, can be es- 
timated by noting that, for a given node j, the maximally 
possible increase in the load before traffic congestion occurs 
is aLj(0). The weight-assignment rule in our traffic proto- 
col, Eq. (0, stipulates that the most probable weight change 
be unity. Now regard a as a control parameter. For a fixed 
amount of change ALj in the load, free flow of traffic is 



a c = ALj/LjiO), 



(3) 



which is independent of the degree variable k [8]. Since the 
load distribution with respect to k is algebraic II 811 . this sug- 
gests that, in order for Eq. (|3} to be meaningful, ALj must 
follow an algebraic scaling law with the same exponent. Since 
the amount of possible weight change is approximately fixed, 
the resulting load change is also fixed. To give an exam- 
ple, we consider a weighted scale-free network of parameters 
N = 3000 and (k) — 4. Initially all nodes are assigned the 
same unit weight. The algebraic load distribution is shown 
in Fig. [3](squares, the upper data set) on a logarithmic scale. 
The algebraic scaling exponent is about 1.5. Next we choose 
nodes of degree k and give them a sudden, unit increase in 
the weight. A recent work shows that for weighted scale-free 
networks, a weight increase of a node typically causes its load 
to decrease |21]. The load change AL as a function of k is 
shown in Fig. [3](circles, the lower data set). We see that on 
the logarithmic scale, AL versus k is parallel to the initial 
load-degree distribution curve, justifying the use of Eq. l|3}. 
Numerically we obtain a* c w 0.37. Since in a realistic sit- 
uation there are more nodes with weights above the uniform 
background value of unity and since the amount of weight 
change can be more than unity, this value of a c is only approx- 
imate. Indeed, direct numerical computations give a c w 0.32. 
The two estimates are nonetheless consistent. An interesting 
observation is that the value of a c is insensitive to the net- 
work parameters like the network size and the degree distri- 
bution, as shown in the two insets in Fig. [3] In particular, for 
N = 1000 and (k) = 4 (inset in the lower-right corner), we 
have a c w 0.39, while for N = 3000 and (k) = 6 (upper-left 
corner), we obtain a c « 0.40. This phenomenon of network- 
parameter independency can be understood by noting that the 
load variation at a node caused by its weight change is mainly 
determined by the probability that optimal paths through this 
node appear or disappear, as a result of the weight change. 
This probability is independent of the network size and the 
degree of the node l21ll . The value of a c , of course, depends 
on the weight-assignment rule and thus on the traffic protocol. 

Can oscillations be expected in realistic networks? To ad- 
dress this question, we test the stability of the Internet at the 
autonomous system level 1 22] . The network comprises 6474 
nodes and 13895 links, the average diameter is (D(0)) w 
4.71, the largest value of the degree is 1460 and the load of 
this node is Lj w 1.97 x 10 8 . By setting a = 0.2, we apply 
perturbation of strength P = 10 x Lj at the largest-degree 
node (to mimic an attack) and let the Internet evolve accord- 
ing to our traffic protocol. The time evolutions of the normal- 
ized diameter (D), of the normalized betweenness (B), and of 
the number of congested nodes nj are shown in Fig. Ela-c). 
Again, persistent oscillations are observed. 

Since network oscillation is a result of the dynamical in- 
terplay between the complex topology and the traffic proto- 
col, this phenomenon is expected to be common for complex 
systems. For simplicity in this paper we adopt the standard 
scale-free random network as the typical model to illustrate 
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FIG. 4: For the Internet at the autonomous system level with capacity 
parameter a = 0.2, evolutions of the normalized diameter (a), of the 
normalized betweenness centrality (b), and of the number of jammed 
nodes. Persistent oscillations of the Internet are observed. 

the oscillation, but extensive numerical evidences have shown 
that this phenomenon is firmly established for complex net- 
works independent of the network parameters like the network 
type, the degree distribution, the average degree, the network 
size, the perturbation size and position, and the specific queu- 
ing scheme of the local node dynamics. We note that although 
the oscillation phenomenon is robust, its time sequence is very 
sensitive to the network details. As a result, the transition from 
periodic oscillations to chaotic oscillations as the capacity de- 
creases is non-smooth. Nonetheless, the trend from regular to 
complex oscillations is still clear. To save space, the details 
will be reported elsewhere l23ll . 

To emphasize the nature of the oscillations, we had avoided 



the fluctuations of capacity parameters and traffic load. While 
these fluctuations will enhance the oscillations, reflected as 
larger oscillation thresholds and larger oscillation amplitudes 
I23B , the deterministic feature of the evolutionary network do 
exemplify the fact that oscillation is an intrinsic property of 
complex systems. This feature makes our model fundamen- 
tally different to the other traffic models in both the temporal- 
and spatial-domain 1 24] . Meanwhile, our model is also dif- 
ferent to the congestion control models where chaotic flux os- 
cillations are observed lo[l4ll . In specific, we consider the 
competitions of simple node dynamics on complex topologies 
while the congestion models investigate the competitions of 
complicated node dynamics on simple topologies. 

In summary, we have discovered that complex network of 
finite capacity can oscillate in the sense that its macroscopic 
quantities exhibit persistent periodic or random oscillations in 
response to external perturbations. While the study is arisen 
as a problem of network security, our findings may have broad 
implications to general traffic networks which, when applying 
the evolutionary weighted model, need further specifications. 
Whereas there can be all sorts of dynamical processes on a 
complex network, our findings suggest that there can be phys- 
ically meaningful situations where the network itself is never 
static but highly dynamic. As a primary model for evolution- 
ary weighted network, the oscillation of macroscopic quanti- 
ties is only one aspect in describing the interplay between the 
topology and the local dynamics, further works in other as- 
pects, such as adaptive networks where the network topology 
updates according to the local dynamics, will be very interest- 
ing. 
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